Transition from plasma- to Kerr-driven laser filamentation 
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While filaments are generally interpreted as a dynamic balance between Kerr focusing and plasma 
defocusing, the role of the higher-order Kerr effect (HOKE) is actively debated as a potentially 
dominant defocusing contribution to filament stabilization. In a pump-probe experiment supported 
by numerical simulations, we demonstrate the transition between two distinct filamentation regimes 
at 800 nm. For long pulses (1.2 ps), the plasma substantially contributes to filamentation, while 
this contribution vanishes for short pulses (70 fs). These results confirm the occurrence, in adequate 
conditions, of filamentation driven by the HOKE rather than by plasma. 
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Filament ation is a self-guided propagation regime typical of high-power lasers, offering spectacular potential 
applications [E[ like rainmaking 0| and lightning control [3. We recently challenged its long-established mechanism 
by measuring the higher order Kerr effect (HOKE) in gases, implying that the non-linear refractive index must be 
written as Anxerr = ^2 n 2jl^, where the non-linear indices n2j are related to the (2j + l) th electric susceptibility 
^/(2j+i) ^he inversion of Anxerr, due to negative 77,4 and ng terms in air and argon, leads to a defocusing Kerr effect 
at an intensity close to those present within filaments d, H|. As a consequence, the HOKE can ensure self-defocusing 
in filaments and balance Kerr self- focusing [Toj], in place of the plasma, especially for short pulses This result 
raised an active controversy flll-[l6| in the lack of direct experimental confirmation. 

Quantitative differences between the predictions of filamentation models including or disregarding the HOKE are 
not sufficient to distinguish between them. The intensity within filaments (~50 TW/cm 2 [17]) is compatible with 
regularization by either the plasma [l8| or the HOKE [To[ , which balance the Kerr self- focusing in the same intensity 
range. Furthermore the electron density is difficult to measure and highly dependent on initial conditions, resulting 
in a wide spread of the reported values from 10 12 to 10 17 cm -3 [18], although the latest measurements range from 
-10 15 cm" 3 [19] to a few 10 16 cm" 3 [ljj. 

In the present Letter, we therefore focus on experimental conditions where qualitatively different behaviors of 
plasma- and HOKE-driven filamentation allows to unambiguously distinguish between them. This discrimination 
proceeds from the intrinsically different temporal dynamics of these non-linear defocusing contributions. While the 
Kerr effect is instantaneous at the timescale of the pulse envelope, the free electrons accumulated throughout the pulse 
duration survive for tens of picoseconds after the laser pulse has passed (20, In an atomic gas like argon, where 
no spatio-temporal modification of the refractive index due to molecular alignment occurs, two pulses separated by 
a few picoseconds can therefore only be coupled if the free electron density left by the first one is sufficient to affect 
the second one. This allows to distinguish between two scenarios. If the plasma is the dominant mechanism for the 
self-guiding of the pump pulse, then the electron density produced in its wake is necessarily affect the probe filament. 
On the contrary, if the HOKE terms are predominant, then the probe filament is insensitive to the presence of the 
pump. 

Based on these different temporal dynamics, we unambiguously observe experimentally the all Kerr-driven fila- 
mentation of ultrashort laser pulses, as well as the transition to a plasma-driven filamentation regime in the case of 
long pulses. This new perspective on the filamentation physics critically revises the discussion of the optimal laser 
parameters for applications ranging from laser-controlled atmospheric experiments to harmonics generation (22I . |23| . 

In our experimental setup (Fig. [1]), two orthogonally polarized co-propagating laser pulses centered at 800 nm, are 
loosely focused by a f=3 m lens in a 2 m-long gas cell filled with argon. We investigate both short pulses (70 fs, 600 /iJ, 
3 bar) and long pulses (1.2 ps, 6mJ, 5 bar), keeping the peak power of each pulse equal to 2.5 critical powers. At this 
power level, each pulse produces a single filament when propagating alone. The probe pulse is temporally delayed 
with regard to the pump by r (r=l ps for short pulses and r=2ps for long pulses). 

The influence of the plasma left by the pump on the probe filament is characterized by observing changes in its 
spectrum. The filament output spectrum is selected in the far-field (z ~2 m) by a pinhole excluding the conical 
emission and most of the photon bath, and analyzed with a spectrometer (Ocean Optics HR4000, providing 0.5 
nm resolution and 14 bits of dynamic range) after separating the two pulses at the cell exit using a Glan-Thomson 
polarizer. To improve the signal-to- noise over the whole considered spectral range, each spectrum is reconstructed 
by assembling data from 3 spectral ranges. For each range, the integration time was adjusted between 2 and 2000 
pulses to ensure the use of the full dynamics of the spectrometer in every spectral region. The resulting spectra were 
then averaged over 20 realizations. The change in the probe spectra induced by the pump pulse is characterized by 
calculating the contrast C(X) = (Si (A) — So (A)) /(Si (A) -fSo(A)) between the spectral densities with (Si) and without 




Figure 1. Experimental setup. Two orthogonally polarized identical pulses with an adjustable delay are focused in an argon- filled 
cell. At the output of the cell, the probe spectrum is selected by a Glan cube polarizer and analyzed with a spectrometer. 
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Figure 2. (a,b) Experimental spectra of the filaments generated by a long (1.2 ps, a) and short (70 fs, b) pulse conveying 2.5 
critical powers, both with and without a pump pulse r=2 ps and r=l ps ahead of it, respectively. This delay ensures that 
the pump pulse can only influence the probe pulse through the defocusing free electrons left behind by the pump. The input 
spectra are also given for reference, (c-d) Contrast between the spectra with and without pump pulse. The probe pulse is 
affected by the plasma left by the pump pulse only in the long-pulse regime. 



(So) the pump pulse at each wavelength A. 

Since the spatial overlap all along the propagation is crucial for the relevance of the measurements, it is optimized 
by maximizing the interference pattern produced by the unfocused pulses both before and after the cell. It is also 
confirmed by the occurrence of multiple filamentation at zero delay, which is set by optimizing frequency doubling 
in a BBO crystal placed before the cell. Moreover, we checked that the alignment is conserved when translating the 
probe pulse from r=0ps to r=lps (resp. r=2ps) by inserting a 200 /im (resp. 400 /im) thick glass window in the 
path of the short (resp. long) pump pulse, delaying it by 1 ps (resp. 2 ps) and checking that the mult i- filament at ion is 
restored. Let us note that the action of the pump pulse on the probe filaments may induce a longitudinal or transverse 
displacement. Such coupling, however, would mean that the plasma strongly affect the filamentation dynamics close 
to the non-linear focus, where a substantial part of the white-light continuum is generated. It would therefore result 
in significant changes in the output spectra. Finally, no broadening is recorded in vacuum, confirming that neither 
input nor output windows of the cell have significant contributions on the spectral broadening. 

As shown in Fig.[2ji, the spectrum of the filament generated by a 1.2 ps long pulse narrows when it is preceded by a 
pump filament r=2ps ahead of it. This coupling demonstrates that the probe pulse propagates through a pre-ionized 
medium with a free electron density providing a significant negative contribution to the refractive index, i.e. a non- 
negligible defocusing term: As expected from the Kerr-plasma balance model, the plasma density generated in the 
filaments substantially contributes to the self-guiding. In this framework, the narrower spectrum of the probe filament 
is easily explained by the supplementary plasma density left by the pump pulse, which decreases the clamping intensity 
ensuring the balance in the probe pulse between Kerr self-focusing and defocusing by the plasma, and consequently 
reduces the efficiency of its spectral broadening. 

Conversely, in the case of short pulses (70 fs, Fig. [2)3), the spectrum of the probe is nearly insensitive to the plasma 
left by the pump. The decoupling between the two pulses separated by a delay (r =lps) too short to allow electron 
recombination unambiguously shows that the plasma density left in the wake of the pump pulse is too weak to signif- 
icantly affect the filamentation process in the considered conditions. Plasma therefore plays no role in filamentation 
of 70 fs pulses, which is instead driven by the HOKE as predicted numerically [Toj]. Note that, if the two pulses 
overlap temporally (r=0), their coupling is restored due to cross-phase modulation, confirming that the two pulses 
indeed overlap both longitudinally and transversely. These results therefore provide the experimental demonstration 
of the transition from plasma-driven filamentation in the case of long pulses to HOKE-driven filamentation for shorter 
pulses, as expected from theoretical considerations 

In order to support this qualitative discussion and provide a closer look on the physical process at play, we numer- 
ically investigate the propagation of the pulses in the present experimental conditions. In a first step, we simulate 
the propagation of the pump pulse. The plasma density left behind this pulse is then used as an initial condition for 
calculating the probe pulse propagation. 

The numerical model considers linearly polarized incident electric fields at a wavelength Ao=800 nm with cylindrical 
symmetry around the propagation axis z. According to the unidirectional propagation pulse equation [24], the scalar 
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Figure 3. (a,b) Spectra of the filaments generated by the long (a) and short (b) pulse, simulated by the full numerical model 
considering the HOKE, (c-d) Contrast between the spectra with and without pump pulse. The probe pulse is affected by the 
plasma left by the pump pulse only in the long-pulse regime. 



envelope e(r, t : z) (defined such that |e(r, z, t)\ 2 = I(r, z, i), I being the intensity) evolves in the frame traveling at the 
pulse velocity according to: 



d z e =i(\J k 2 (uj) - k\ — k'uj)e 



iu 2 ~ uj ~\ W 

-Pnl — - — -J 



where c is the velocity of light in vacuum, uo is the angular frequency, /c(cj)=n(cj)cj/c, k' its derivative at cjo=27tc/Ao, 
n(oS) is the linear refractive index at the frequency lj, k± is the spatial angular frequency. Pnl is the nonlinear 
polarization, J is the free-charge induced current and a is the nonlinear losses induced by photo-ionization. / denotes 
simultaneous temporal Fourier and spatial Hankel transforms of function /: / = J °° J^° oo rJo(k±r)f(r^t)e %UJt dtdr J 
where Jo is the zeroth order Bessel function and / = £,Pnl, or a. Pnl is evaluated as Pnl = ^2 n 2j\^\ 2 ^^ 
where ri2j are the j th -order nonlinear refractive indices as measured in The current is evaluated as J = ^-(^ e + 
iuj)ps/(v 2 -fa; 2 ), where e and m e are the electron charge and mass respectively, eo is the vacuum permittivity, v e is 
the effective electronic collisional frequency, and p is the electron density. Finally, a = W(\e\ 2 )Ui (p a t — p) /(2|^| 2 ), 
p a t is the neutral atoms density, iy(|e| 2 ) is the photoionization probability modeled by the PPT (Perelomov, Popov, 
Terent'ev) formulation, with ionization potential Ui. 

The propagation dynamics of the electric field is coupled with the electron density p, calculated as (H 

d t p = W{\e\ 2 ) (p at -p) + ^p\s\ 2 - (3p 2 , (2) 

where f3 is the electron recombination rate and a is the inverse Bremsstrahlung cross-section of argon, also accounting 
for avalanche ionization. The output spectrum is integrated over 2 mrad around the beam center to match the 
experimental conditions. 

The full model including the contribution of the HOKE together with the plasma response reproduces remarkably 
well the experimentally observed behavior. In particular, as displayed in Fig. [3j the spectral broadening of the long 
probe pulse is reduced when the pump pulse precedes it, in a ratio comparable with the experimental data, while the 
filament generated by the 70 fs probe pulse is unaffected by the pump pulse. In contrast, disregarding the HOKE 
(Fig. [4]) would lead us to expect that the filament output spectrum generated by probe pulses of any duration should 
be affected by the presence of the pump pulse, in contradiction with the experimental results. Furthermore, this 
truncated model inadequately predicts the wings of the spectral broadening, even for a single pulse. Comparing 
both Figs. [3] and 2] with Fig. [2] clearly illustrates the need to consider the HOKE in numerical simulations of laser 
filamentation, even in the plasma-driven, long-pulse filamentation regime. 

This need is illustrated by comparing the electron densities predicted by both models. While the truncated model 
yields 10 17 cm -3 for both 1.2 ps and 70 fs pulses, the full model yields 4xl0 16 cm -3 and 10 15 cm -3 , respectively. 
This strong dependence of the plasma contribution on the pulse duration explains both the contrasted behaviors 
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Figure 4. (a,b) Spectra of the filaments generated by the long (a) and short (b) pulse, simulated by the truncated numerical 
model disregarding the HOKE, (c-d) Contrast between the spectra with and without pump pulse. The probe pulse is affected 
by the plasma left by the pump pulse whatever the pulse duration. 



observed in our experiment between the plasma- and HOKE-driven filamentation regimes, but also the wide spread 
of the experimentally measured electron densities in filaments [l3[ fl9| . Furthermore, we can estimate the relative 
contributions of HOKE and plasma to defocusing by considering the ratio £ = AnHOKE/An p i asma of the HOKE- to 
plasma- induced refractive index change. For short pulses (70 fs), this ratio keeps well above 1 all along the filament 
length (£ > 39), illustrating the negligible contribution of plasma to the filamentation process. In contrast, for 1.2 ps 
this ratio keeps close to 1 (£ > 0.72, except for a spike with £ = 0.24 at the very filament onset), confirming that, 
while plasma provides the major defocusing contribution, the HOKE are far from negligible even in these conditions. 

As a conclusion, a pump-probe experiment allowed us to unambiguously observe experimentally the theoretically 
predicted HOKE-driven filamentation for ultrashort pulse sflG | , as well as the transition from this regime to the long- 
known plasma-driven filamentation regime for long pulses [9] . This transition is similar to that observed in the context 
of high-harmonic generation (HHG), where the use of too long pulses results in gas ionization instead of HHG (23j . 
Furthermore, comparing our results with numerical simulations show that, even in the plasma-driven filamentation 
regime of the present work, the contribution of the HOKE to the propagation dynamics cannot be neglected. This 
finding provides a better understanding of filamentation, and therefore allows to improve its modelling. It further 
confirms the relevance of the measured HOKE |8, 9], with implications ranging from spectral broade ning in optical 
fibers [25] to the generation of few-cycle pulses jUj], atmospheric applications |H|43|? or fermionic light (27f . 
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